// Коэффициент отражения принят равным единице, фаза при отражении %pi.
function [Ez, Z2]=ant_ran_mp_ground(Ez,Z2,Na,L,lambda,H,EA, td,ha)  //
exec('Rn1.sci');
R=Rn1(EA); 
[u,v]=size(EA);
Tau=3.3*ha*abs(sin(EA(:,1))); // № строки - № точки отражения,  № элемента в строке:
                                                                                // временной сдвиг для сигнала с этим номером

Tau=round(1e-9*Tau/td);
//rem = Tau*1.602+1/2 - fix((Tau*1.602+1/2)./1).*1;
r=rem(Tau*1.602+1/2,1);
ph=exp(%i*r*2*%pi);//
if Na==3
    exec('antset_quad.sci');
    Pa=antset_quad(Na,L,H);      // Орты базовых линий; Na=3;  (строка - орт базы) 
end;
if Na==6
    exec('antset1.sci');
    Pa=antset1(Na,L,H);                // Орты базовых линий; Na=6;  (строка - орт базы) 
end;
if Na==1,
    Pa=[L,0,0];
end;
PaR=Pa*R';                           // матрица проекций ортов  линии распространения сигналов на ант. базы.
                                                // № строки - это № антенны, № элемента в строке -- № сигнала
PH=exp((2*%pi*%i/lambda)*PaR);            // матрица фазовых сдвигов;  № строки - это № антенны, 
                                                                     //№ элемента в строке -- № сигнала
exec('circshift.sci');
                                                     // i - № точки отражения,
for  k=1:u                                                 //  k --№ сигнала
//     A = Ez(k,:)';
//     shift = Tau(k);
//     cs = [A(&,:);A(shift:$-1,:)];
//     Ez1(k,:)=0.7*cs'*ph(k);  //Kmp(1)*
    Ez1(k,:)=0.7*circshift(Ez(k,:)',0,Tau(k))'*ph(k);  //Kmp(1)*
    Ez(k,:)=Ez(k,:)+Ez1(k,:);                     
    for    j=1:Na,                                         // j --№ антенны
        Z2(j,:)=Z2(j,:)+Ez1(k,:)*PH(j,k);          
    end;
end;

endfunction





